Electric field in stationary superconductors 
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It is generally accepted that vortex core is charged, what illustrates that even in stationary super- 
conductors electric field may be present. Vortex structure and other properties of superconductors 
are usually calculated in the framework of the Ginzburg-Landau theory, which does not cover elec- 
tric field. We show that the generalization of the GL theory due to Bardeen allows one to derive a 
third GL equation for the electrostatic potential. Since the Bardeen theory applies to all tempera- 
tures, the presented theory enables one to calculate charge profiles of the vortex under quite general 
conditions. The theory is consistent with the BCS-Gorkov results. 



I. INTRODUCTION 



II. GORTER-CASIMIR TWO FLUID MODEL 
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As was pointed out by Adkins and Waldram Q and 
later followed by others, see e.g. Q, a difference of the 
chemical potential in a superconducting versus normal 
state, is compensated by the electrostatic potential. If 
the gap is modulated in space, say due to supercurrents, 
this potential results in an internal electric field. This 
idea has been used later to estimate the charge of vortices 
HQ. For superconductor of the first kind, the presence 
of nonzero electric field in the superconductor carrying 
tranport current was first experimentally proved by Bok 
and Klein [Q who measured the surface charge. For re- 
cent measurement with this method see e.g. A new 
experiment, which allows one to measure the electric field 
in the bulk of the superconductor of the second kind in 
the mixed state, has been performed by Kumagai, Noza- 
kii and Matsuda 0. 

A vortex structure is usually calculated in the frame- 
work of the Ginzburg-Landau (GL) theory, which is ap- 
plicable only for temperatures close to the critical tem- 
perature T c and in its standard formulation does not 
cover the electric field. The electric field can be evalu- 
ated, however, from the Poisson equation, called by Jake- 
man and Pike the third GL equation, to reflect that 
the source term (density of unscreened charge) is a func- 
tional of the GL function. 

Various expressions for the source term of the third GL 
equation can be found in the literature Q-QB-M, most 
of contributions are of similar amplitude and seem to be 
merely alternative approximations of the same mecha- 
nism. This is not, however, the case. There are at least 
three distinguishable mechanisms which create the elec- 
tric field in superconductors. To clarify this point, we 
derive the third GL equation with all these mechanisms 
from a simple phcnomenologic theory of the GL type. 
Not to be restricted to temperatures close to T c , we adopt 
Bardeen's extension of the GL theory |H] based on the 
Gorter-Casimir model |Ol. 



The free energy of the normal state metal without elec- 
tric and magnetic fields reads 



?n = U - \lT 2 



(1) 



where U is internal energy including lattice vibrations 
and 7T is the electronic specific heat. According to the 
Gorter-Casimir two fluid model, the free energy of a su- 
perconducting state without fields can be written as 
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T S = U- -7T c a w - - 7 T 2 v / T _ ^. 
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Equilibrium value of the order parameter u> makes J- s 



minimum. From the condition = it follows that 



= 1-t 4 



(3) 



where t = T/T c is the reduced temperature. At T = the 
equilibrium value oj cq — 1 and T n —!F$ = \ jT^ = \[iR 2 , 



where H c is the thermodynamic critical field. At T = T c 
the equilibrium value of lu is zero and J- s = T u = 
U — \^fT 2 as it should. According to its temperature de- 
pendency, the order parameter ui can be identified with 
the square of the GL effective wave function -0, normal- 
ized to the superfluid fraction, |-0| 2 = where n s is 
the superconducting charge carriers density and n is the 
total density of charge carriers. 

Near the critical temperature the order parameter u> is 
small. Using the expansion a/1 — w ~ 1 — \^ — \ uj2 the 
Gizburg-Landau free energy, a |i/'| 2 -|-i/3 j^l is recovered, 
where the parameters a and (3 in the used normalization 
are a = —2^,H 2 (l — t)/n and — fiH 2 /n 2 . Accordingly, 
one can see the free energy of Ginzburg and Landau as 
an asymptotic form of the more general Gorter-Casimir 
model. 



III. BARDEEN'S EXTENSION OF THE GL 
THEORY 



To extend the region of applicability of the Ginzburg 
Landau theory, Bardeen (Til] replaced the Ginzburg 
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Landau polynomial free energy by the free energy due 
to Gorter and Casimir. The free energy then reads 
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We note that the last two terms, which describe the elec- 
trostatic field, have not been assumed by Bardeen [jy] . 
The magnetic and electric fields are given by standard 
definitions, B = VxA and E = —Vip. 

The state of the system is given by the minimum of 
the total free energy. Accordingly, variations of T with 
respect to the vector potential A, the scalar potential ip, 
the GL function ip and the electron density n have to 
vanish. 

Variation ~ — leads to extended Ginzburg-Landau 

dtp ° 

equation, sec Q 
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Variation |^ = leads to the Maxwell equation, 



V x B 



-lien — Ret/jfiKV + e*A)ip. 
m* 



(6) 



In principle, density n of the total charge is modified 
by internal electric fields. Accordingly, these equations 
should be treated together with the Poisson equation pre- 
sented below. Since both equations depend on the total 
density n, not on the deviation from the homogeneous 
value, it is clear that for calculating the magnetic prop- 
erties one can safely neglect the effect of the small electric 
field on the density. Within this approximation, the set 
(§) and (§) is closed. 
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E F (n) « E% + Sn 

ployed the Poisson equation. Accordingly, equation (0) 
can be viewed as the Poisson equation with screening. 

The rest of terms in the right hand side of (Q) are source 
terms of the third GL equation corresponding to various 
mechanisms. The second term is the internal pressure 
due to the dependence of the condensation energy on 
the total density n. The third term is thermoelectric 
field of the normal metal reduced by factor y/l — |^| 2 . 
The fourth term of is the non-local Bernoulli potential 
fLOf but reduced by \ip\ 2 in the spirit of the so called 
quasiparticle screening |l3[ . The correction due to the 
density dependency of the Cooperon mass, cx d ^[™ H , has 
been derived by Rickayzen II but it is ussually omitted. 



V. COMPARISON WITH THE BCS RESULTS 



Using the BCS formula for the condensation energy 
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£con = \lT 2 c = \liHl = l -NA 2 , 



(8) 



with Ao = 2Uu)d exp(— 1/A/V), and the relation of the 
gap to the GL function, A = AoV', one recovers from the 
second term of (Q) the BCS expression, 
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The first expression has been derived e.g. in The 
second form uses the parabolic band approximation, 
^ = _i , Together with the estimate, AfV ~ 1, it 

On It^F ° r—i 

gives the formula used e.g. in M. 



IV. THIRD GL EQUATION 



VI. CONCLUSION 



Variation = leads to the Poisson equation, 

yields the third 
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—eV <p — p. Finally, variation -j— 
GL equation, 
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The first term of (|7|) represents the screening on the 
Thomas- Fermi lenght ||], A^ F = 2 Afe 2 > w here JV is 
the single-spin density of states. In deriving this term 
we have used the linear response approximation, = 



Starting from the combination of the Gorter-Casimir 
and Ginzburg-Landau free energies further extended by 
the energy of the electrostatic field, we have derived a 
set of three GL equations. This set includes the ex- 
tended GL equation (|5|), the Maxwell equation (||) for 
the magnetic field and the Poisson equation (^) for the 
electrostatic potential. In particular limits, the theory 
reproduces previous results. 
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